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Abstract 

Two paradoxical aspects of the prevailing kinetic equations are pre- 
sented. One is related to the usual understanding of distribution function 
and the other to the usual understanding of the phase space. With help of 
simple counterexamples and direct analyses, involved paradoxes manifest 
themselves. 

PACS number: 51.10.+V. 



When the Boltzmann equation, thought of as the first equation of the prevail- 
ing kinetic theory, came out, strong doubts arose, of which many were related 
to the fact that whereas Newton's equations themselves were time-reversible 
the Boltzmann equation, as a consequence of Newton's equations, was time- 
irreversible [1]. As time went by, particularly after the "rigorous" BBGKY the- 
ory was formulated in the middle of the last century [2— 8], the philosophical 
concern of the time-reversal paradox gradually faded out, and it was believed 
that the ultimate understanding of the related issues had been completely es- 
tablished, at least in the regime of classical mechanics. 

However, relatively recent developments of mathematics and physics seem 
to have brought new elements into the picture. In particular, the studies of 
fractals[9, 10] reveal that there can, at least in the mathematical sense, exist 
functions that have structures of self-similarity even at the infinitesimal level. 
Such functions are intrinsically discontinuous and cannot be described by usual 
differential apparatuses. Along this line, an increasing number of scientists are 
surmising that if any similar structures are found in realistic gases, some of 
conventional concepts in the standard theory need to be revised significantly. 

In connection with this, we wish to present and discuss two relatively un- 
known aspects of the standard framework of kinetic theory. Firstly, it will be 
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shown that a realistic gas contains, almost always, a significant amount of par- 
ticles whose distribution function does not keep invariant in the six-dimensional 
phase space and cannot be regarded as a continuous one. Secondly, it will be 
illustrated that there are inherent difficulties in formulating the particles that 
enter and leave an infinitesimal phase volume element during an infinitesimal 
time, which suggests in the sense that the phase space is more sophisticated 
than the customary thought assumes. By removing abstractness of the matters, 
all related paradoxes become surprisingly simple and straightforward. We are 
now convinced that if those paradoxes had been unveiled at the very berginning, 
kinetic theory would have been renewed several times. 

The basic core of the classical kinetic theory says that if particles of a gas do 
not interact with each other, the distribution function describing them satisfies 
the collisionlcss Boltzmann equation 

df df F df n 

at or m av 

which is, according to well-known textbooks[ll, 12], equivalent to the path- 
invariance of collisionless distribution function 



dt 



= 0, (2) 

r(t),v(t) 



where the subindexes r(t) and v(f) imply that that the differentiation is taken 
along a particle's path. Conceptually speaking, equation (2) can be interpreted 
as saying that such a gas is "incompressible" in the phase space. The standard 
theory further states that, if collisions between particles are more than negli- 
gible, a certain type of collisional operator needs to be introduced to replace 
the zero term in (1). After having the collisional operator due to Boltzmann 
himself, equation (1) becomes 

% + + I n lf(v')f(vi)-fMf(vi)}ua(n)dndv u (3) 

where the first term on the right side represents particles entering the unit phase 
volume element and the second term particles leaving the unit phase volume 
element. 

Although the formulation briefed in the last paragraph seems stringent and 
has been accepted unanimously in the community, it is not truly sound. Ironi- 
cally enough, even a simple glance at the form of the Boltzmann equation (3) 
offers intriguing things to ponder. On the left side there is a symmetry between 
r and v in terms of differentiation operation; whereas on the right side the po- 
sition vector r serves as an inactive "parameter" and all the integral operations 
are performed in the velocity space. (The solid angle is defined in terms of 
the velocities.) This disparity, while seeming a bit too schematic to be fully 
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convincing, should serve a motivation for us to investigate the whole issue more 
carefully and more thoroughly. 

We first look at whether or not the collisionless Boltzmann equation (1) and 
the path-invariance theorem (2) make sense as they intend to. 

In the mathematical sense, the picture provided by the path-invariance is 
quite clear and simple. Consider a sufficiently small volume element 

AxAyAzAv x Av y Av z (4) 

moving together with a certain particle in the phase space; the theorem asserts 
that the particle density within the element (in number herein) does not increase 
or decrease. In well-known textbooks, the theorem is "proved" by applying the 
rigorous Jacobian approach[ll — 13]. Here, for the purpose of this paper, we 
wish to illustrate the theorem and its derivation in an intuitive way. Fig. 1 de- 
scribes what takes place in the two-dimensional x-v x subspace. As time passes 
from t = to t — T, the particles distributed inside the rectangle AxAv x in 
Fig. la will be distributed inside the parallelogram in Fig. lb. By excluding all 
external forces (just for simplicity), it is obvious that the area of the parallel- 
ogram is equal to that of the rectangle and the average particle density within 
the moving volume element keeps invariant. On the understanding that parti- 
cles distribute continuously within AxAv x and the size of AxAv x can shrink to 
zero, the invariance of the average particle density can certainly be interpreted 
as the path-invariance of distribution function. The discussion above, though 
simplified somewhat, reflects the essence of the full theorem accurately. 

It seems, at this stage, that the validity of the path-invariance (2), together 
with the picture of incompressible fluid in the phase space, is so solid and so clear 
that we should discuss it no more. But, somehow, the truth is not that simple: 
there exist many cases in that gases behave very differently. To get an immediate 
idea about such behavior, let's look at the special arrangement shown in Fig. 2, 
where particles with definite velocity v strike a convex solid surface and then 
get reflected from it elastically. (Later on, more general models about collisions 
between particles and boundaries will be adopted and examined.) Following a 
moving particle and counting particles in a definite position range dr = dxdydz 
and in a definite velocity range dv = dv x dv y dv z , we find out a clear-cut fact 
that the reflected particles simply obey 



dt 
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That is to say, these particles diverge in the phase space. In a similar way, we 
can readily see that if the incident particles fall upon a concave solid surface, 
the reflected particles will, in some region, obey 



dt 



> 0. (6) 
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This tells us that the particles converge. 

A sharp question arises. How can the path-invariance theorem, considered 
as the very core of all the kinetic equations, possibly suffer from such simple 
and direct counterexamples? By reviewing the derivation of the Boltzmann 
equation, we realize that equations (1) and (2) hold only for perfectly continuous 
distribution functions, while the diverging and converging particles presented 
above are related to none of them. 

To see the point more vividly, let's go back to the two-dimensional phase 
space, namely the x-v x space, and investigate the evolution of the particles 
marked with the dotted curve in Fig. 2. In the schematic sense, we may say 
that at t = all these particles are distributed along the diagonal of AxAv x 
shown in Fig. 3a, rather than in AxAv x uniformly. After a while, at t = T, this 
diagonal is stretched and becomes longer as shown in Fig. 3b. If we set up a 
small, but definite, phase volume element and let it move together with one of 
the particles, we will certainly find that the particle density within it decreases. 
(For discontinuous distributions, if we allow the shape and size of measuring 
volume element to vary, we can get any value, drastically from zero to infinity, 
for the particle density.) 

Now, it is in order to examine how a realistic boundary surface S "reflects" 
realistic particles that come from all directions and with different speeds. Exper- 
imental facts inform us that the reflection cannot be truly elastic and stochastic 
and dissipativc effects must play a certain role [14]. To express particles "pro- 
duced" by such surface, it is proper (particularly if the gas is a rarefied one) to 
define the instant emission rate a on an infinitesimal surface dS in such a way 
that 

adtdSdvdfl (7) 

represents the number of particles ejected by the surface element in the speed 
range dv and in the solid angle range dQ during the time interval dt. By al- 
lowing the emission rate a to be described in probability and to have certain 
dependence on velocities of incident particles, it can be said that the stochastic 
and dissipative nature of the reflection has been included. Then, divide S into 
N surface elements. Referring to Fig. 4, we find that for the ith surface element 
(AS)i the reflected particles are like ones emitted from a "point particle source" 
at roi, the position vector of (A5)», and the relevant distribution function at a 
point r in the reflection region takes the form 

f(r, v, Q) = a{AS ^ 3 Vi(v)5(Sl n T -r 0i ), (8) 
|r - rojfir 

where Vi{v) is a certain function of v and 5(Q — O r _ roi ) is the S- function defined 
on the solid angle in the velocity space. It is very obvious that, regardless of 
the forms of a and V(v), this distribution function diverges in the phase space. 
Another interesting point about the distribution function is that it is perfectly 
continuous in terms of r and v and it is like a function defined on a single point in 
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terms of tt. Based on this observation, we are tempted to say that the function 
is on a variable domain of 4 + e dimensions. The total distribution function at 
r associated with all reflected particles from S is 



N 



f(T,V,Sl)=Y, 



|r - r i| 2 u 3 



Vi(v)6(n-ci r . roi ). 



(9) 
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Though the distribution function expressed by (9) can, as N — > oo, be regarded 
as a sum of an infinitely large number of (4 + e)— dimensional functions, it is 
not a continuous function defined in the phase space. Actually, with help of this 
expression, we can analytically or numerically prove that the particle number 
within a given moving volume element drdv is capable of decreasing, increasing 
or keeping invariant. As a limiting case, we may assume S to be relatively small, 
or the distance r — r i | to be relatively large, and find that expression (9) is 
reduced to expression (8) and the particles related to it always diverge in the 
phase space. 

The above investigation, though formally simple, provides us with a very 
new picture on gas dynamics. Rather than as a continuous medium or an in- 
compressible fluid in the phase space, a gas should be considered as a special 
collection of discrete particles, whose distribution function can change from con- 
tinuous one to discontinuous one, as well as from discontinuous one to continu- 
ous one. Another notable fact, which may very much interest scientists carrying 
out practical studies, is that the changeovers aforementioned occur dramatically 
near interfaces between fluids and solid boundaries. 

We now turn our attention to the validity of the collisional operator on 
the right side of (3). A widely accepted concept is that in order to formulate 
collisional effects we are supposed to focus ourselves on a fixed six-dimensional 
phase volume element and study how particles leave and enter the element due 
to collisions. In what follows, it will be shown that such concept is, much to our 
surprise, deceptive. 

Firstly, particles leaving a phase volume clement drdv during a time interval 
dt are of interest. According to the standard theory, if a particle collides with 
another within the volume clement drdv during dt, it should be considered as 
one that leaves the volume clement during the time interval because of the ve- 
locity change caused by the collision (see Fig. 5). This intuitive, seemingly very 
reasonable, picture can be challenged in the following way. In deriving kinetic 
equations, a necessary step is to let dt, dr and dv approach zero independently. 
If it is assumed, in the limiting processes, that the length scale |dr| is much 
smaller than |vdi|, then virtually all the particles, initially within drdv, will 
leave drdv at the end of dt, irrespective of suffering collisions or not. If we still 
want to say that the standard consideration, concerning how many particles 
stay inside drdv without involving collisions and how many particles get out of 
drdv with collisions involved, holds its significance, we have no choice but to 
assume that |dr| >> |vdi|. An unfortunate fact is that no sound reason can be 
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found out for that we can prefer this assumption to its converse. 

Secondly, particles entering a phase volume element drdv during dt are under 
examination. In the standard treatment, two beams with velocities v' and 
v[ are assumed to collide within a volume element dr during dt and to give 
contribution to the particles expressed by f(r,v,t)drdv, as shown in Fig. 6. 
Ironically enough, this treatment involves not one but many paradoxes. As 
one thing, particles produced by collisions in a small spatial region, like ones 
emitted from a point particle source, will diverge in the phase space and cannot 
be treated as an ordinary contribution to f(t,r,v)drdv. As another thing, we 
again need to let dt, dr and dv approach zero. If |dr| << |vdi|, all produced 
particles will "instantly" leave drdv and only an insignificant fraction of them 
can be treated as a contribution to f(r,v,t)drdv. See Ref. 16 to get more 
paradoxes. 

In summary, two paradoxical aspects of the standard kinetic theory have 
been presented. The first aspect is related to the distribution function. A 
tacit assumption of the existing kinetic equations is that distribution functions, 
though describing discrete particles, must be mentally and practically contin- 
uous. The discussion of this paper, however, shows that distribution functions 
of realistic gases have, in general, complex local structures and they cannot be 
described by continuous distribution functions. The second aspect is related to 
the phase space. A usual concept in the standard theory is that the position 
space and the velocity space can be separated mentally and practically. The 
discussion of this paper, however, shows that whenever we investigate the time 
dynamics of particles in the velocity space we should keep an eye on what takes 
place in the position space, and vice versa. 

A variety of fundamental questions can be raised, of which many are beyond 
the scope of this brief paper. In some of our recent works, we make more analyses 
and put forward alternative approaches[15, 16]. With help of a development 
in quantum mechanics [17], some of the discussion in this paper can also be 
extended to the regime of quantum statistical physics. 

This paper is supported by School of Science, BUAA, PRC and by Education 
Ministry, PRC. 
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Fig. 1: A moving phase volume clement (a) at an initial time t = 
and (b) at a later time t = T. 
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Fig. 2: Schematic of reflected particles from a solid surface. It is 
interesting to note that whether or not the involved collisions are 
elastic will not alter this picture significantly. 
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Fig. 3: Schematic of how the particles marked in Fig. 2 spread 
in the x — v x space, (a) These particles are distributed along one 
diagonal of the rectangle AxAv x at t = 0, and (b) the diagonal is 
stretched at t = T. 



8 




Fig. 4: Particles reflected from a small surface element of realistic 
boundary. 




Fig. 5: A particle involving a collision in a phase volume element. 




Fig. 6: Two beams of particles collide with each other and produce 
particles with velocity v. 
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